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ABSTRACT In telomerase-deficient yeast cells, like equivalent mammalian cells, telomeres shorten over 
many generations until a period of senescence/crisis is reached. After this, a small fraction of cells can 
escape senescence, principally using recombination-dependent mechanisms. To investigate the pathways 
that affect entry into and recovery from telomere-driven senescence, we combined a gene deletion dis- 
rupting telomerase [estIA) with the systematic yeast deletion collection and measured senescence charac- 
teristics in high-throughput assays. As expected, the vast majority of gene deletions showed no strong 
effects on entry into/exit from senescence. However, around 200 gene deletions behaving similarly to 
a rad52A estIA archetype {rad52A affects homologous recombination) accelerated entry into senescence, 
and such cells often could not recover growth. A smaller number of strains similar to a riflA estIA archetype 
{hf1A affects proteins that bind telomeres) accelerated entry into senescence but also accelerated recovery 
from senescence. Our genome-wide analysis identifies genes that affect entry into and/or exit from telomere- 
initiated senescence and will be of interest to those studying telomere biology, replicative senescence, cancer, 
and ageing. Our dataset is complementary to other high-throughput studies relevant to telomere biology, 
genetic stability, and DNA damage responses. 
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In most eukaryotes, including mammals and yeasts, telomerase helps 
maintain telomeric DNA and solve the chromosome end- replication 
problem (Gomes et al. 2010). Telomeric DNA recruits a number of 
different telomere-binding proteins to make functional telomeres. 
Dysfunctional telomeres can stimulate DNA damage responses, lead- 
ing to genetic instability, cell senescence, or apoptosis (Jain and Cooper 
2010; Lydall 2009; Stewart and Weinberg 2006). 
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In most human somatic cells, there is insufficient telomerase to 
maintain telomere length indefinitely through cell divisions, and 
consequently, telomeres shorten and ultimately become nonfunctional. 
Telomere shortening is thought to be a double-edged sword in meta- 
zoan cells (Gomes et al. 2010; Stewart and Weinberg 2006). On one 
hand, the limited lifespan of functional telomeres limits cell divisions and 
therefore protects against uncontrolled cell division in cancer (Artandi 
and Depinho 2009; Lin et al. 2010). On the other hand, limiting cell di- 
vision prevents tissue renewal and repair and thus contributes to ageing. 

Telomere structure and function need to be maintained in cells 
with unlimited proliferative potential. In the majority of cancers, 
structure and function are preserved via the expression of high levels 
of telomerase, permitting unlimited and uncontrolled cell division. 
However, in around 15% of cancers, telomeres are maintained by al- 
ternative lengthening of telomeres (ALT) mechanisms (Henson et al. 
2002). Therefore a fuller understanding of the pathways allowing eukary- 
otic cells to tolerate dysfunctional telomeres or to maintain functional 
telomeres, either in the presence or absence of telomerase, is critical 
for understanding processes underlying carcinogenesis and aging. 
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In budding yeast, telomerase is constitutively expressed and telo- 
mere length is stable. However, telomerase can be deleted, giving yeast 
cells a life-cycle much like mammalian somatic cells; their telomeres 
shorten, they enter into a period of senescence/crisis, and at low rates, 
the cells recover from senescence using telomerase-independent mech- 
anisms to maintain telomere length and function (Lundblad and 
Blackburn 1993). Telomerase-deficient yeast cells are therefore an ex- 
cellent model system to study the pathways that allow eukaryotic cells 
to tolerate dysfunctional telomeres or to maintain telomeres in the 
absence of telomerase, and are thus useful for studying pathways 
relevant to human carcinogenesis. Telomerase-deficient yeast cells 
typically survive by using recombination-dependent mechanisms to 
maintain telomeres and are a model for ALT pathways in human cells 
(Lendvay et at 1996; Lundblad and Blackburn 1993). However, re- 
combination-independent mechanisms for chromosome end protec- 
tion also occur in budding and fission yeasts (Jain et at 2010; Lee et at 
2008; Maringele and Lydall 2004b). 

Although a number of yeast genes, principally genes affecting 
homologous recombination, have been identified that affect entry/exit 
from telomere driven senescence (Abdallah et at 2009; Azam et at 
2006; Chang et at 2011; Enomoto et at 2002; Enomoto et at 2004; 
Ijpma and Greider 2003; Le et at 1999; Lee et at 2007; Lydeard et at 
2007; Maringele and Lydall 2004a; Mceachern and Haber 2006), so far 
no genome-wide analysis of these processes has been reported. Here 
we describe the first genome- wide survey to systematically determine 
the effects of budding yeast gene deletions on the growth of cells with 
telomere defects. Through these systematic approaches in budding 
yeast, we identify many yeast genes that affect telomere- driven senes- 
cence, orthologs of which may play important roles during carcino- 
genesis or aging. 

MATERIALS AND METHODS 

Growth media, synthetic genetic array, and strains 

An estlr.natMX strain (DLY5026) was combined with the genome- 
wide knockout collection using SGA as previously described (Addinall 
et at 2011). W303 genetic background strains were cultured in yeast 
extract, peptone, dextrose, and adenine (YEPD) media. 

Culture passage 

Solid passage: 1536 colony solid agar-to-solid agar pinning was 
performed on a Biomatrix BM3-SC robot (S&P Robotics Inc., Tor- 
onto, Canada) using an 0.8 mm diameter pin tool. 

Liquid passage: Selected estlr.natMX yfg::G418 double mutants were 
cherry-picked after SGA. Inoculation from solid to liquid media was 
performed on the Biomatrix BM3-SC robot using a 96-pin (1 mm 
diameter) pintool. 96-well liquid culture plates were incubated at 30°C 
for 22 hr without shaking. A Biomek FX robot (Beckman Coulter 
Limited, High Wycombe, UK) was used to vortex the cultures to 
resuspend and then inoculated into new sets of 96-well plates con- 
taining 200 jjlI media in each well. A 96-pin (2 mm diameter) pin tool 
(V&P Scientific, Inc., San Diego, CA, USA) was used, transferring 
approximately 3 ul of cell culture at each passage. In a similar pro- 
cedure, resuspended cell cultures were diluted in 200 ul sterile H 2 0 
and vortexed before the diluted cultures were spotted onto solid agar 
plates. Agar plates were incubated for 48 hr at 30°C before being 
photographed. 

Culture size estimation: Solid agar plates were photographed on an 
splmager (S&P Robotics Inc.) and analyzed using Colonyzer (Lawless 



et at 2010). Culture size variances were stabilized using a Box-Cox 
transformation with a power parameter of 0.6 as implemented in the 
R statistical package (R Development Core Team 2011). 

Mean density profile similarity: The similarity of the MDP shape of 
each strain to that of three example archetypes was estimated using 
Pearson's correlation coefficient (r) between strain MDPs and the 
MDP of each archetype, using the cor function in the R statistical 
package (R Development Core Team 2011). As well as calculating r, 
absolute differences in MDP were estimated using the root mean 
square (RMS) difference between profiles. The area under each 
MDP (representing total strain fitness across all passages) was also 
estimated as the sum of mean culture size across all 22 passages. Lists 
ranking the MDP for each strain according to its correlation coeffi- 
cient were generated, with the strains having MDP shape most similar 
to the archetype (and the archetype itself) at the top of the list (Figure 
4, and supporting information, File S5, File S6, File S7, and File S8). 
Note that the absolute fitness of a comparator strain does not have to 
be similar to that of the archetype for their MDPs to correlate well. 
Two strains, one sick and one healthy, can have similarly shaped 
MDPs that correlate well. It is for this reason that some of the heat 
maps in Figure 4 are perceptibly stripy. 

Fraction of recovered repeats: Recovered fraction for each genotype 
was estimated by counting the number of repeated cultures for each 
strain and calculating the fraction of those repeats which had recovered 
to a size of 1000 or greater by passage 19 (typically in the range from 0/8 
to 8/8). For classes similar to each of the three archetypes, the average 
recovered fraction of class members was calculated. 

Telomere length detection 

Southern hybridization to examine telomere length and structure was 
performed similarly to that previously described (Maringele and 
Lydall 2004a). Genomic DNA was extracted, digested with Xhol, then 
run overnight on a 1% agarose gel at 1 V/cm. Southern transfer 
and detection was performed using DIG-High Prime Labeling and 
Detection Kit (Roche) and visualized on a FUJI LAS4000. Telo- 
meric Y'-TG probes (pHT128; Tsubouchi and Ogawa 2000) were 
labelled using DIG-High Prime DNA Labeling and Detction Starter 
Kit II. 

RESULTS 

Measuring senescence by high-throughput, 
multipassage assays 

To identify gene products that affect how telomerase-deficient cells 
enter into and recover from senescence, we used the synthetic genetic 
array (SGA) method to combine a systematic genome-wide collection 
of gene deletions with an estlA mutation defective in telomerase 
(Costanzo et at 2010). As recovery from senescence is a rare and 
stochastic process (Lundblad and Blackburn 1993; Nautiyal et at 
2002), we compared two high-throughput methods of measuring se- 
nescence. These two methods were analogous to low-throughput liq- 
uid culture or serial tooth-picking approaches (Lundblad and 
Blackburn 1993; Maringele and Lydall 2004a). In a liquid procedure 
we passaged yeast strains in 96-well plates, while in a solid agar pro- 
cedure we passaged 1536 colonies per plate by pinning directly onto 
agar (Figure 1). The liquid procedure allowed greater dilution and 
therefore a greater range of divisions at each passage, but throughput 
was less efficient than in the pinning procedure. The liquid procedure 
had a higher risk of culture contamination. 
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Figure 1 Two methods for ge- 
nome-wide analysis of entry 
into and recovery from senes- 
cence. An est7::natMX strain 
(DLY5026) was crossed to the 
deletion mutant collection and 
double mutants were selected 
using SGA procedures. (A) Dou- 
ble mutants were passaged in 
96-well plate liquid cultures, 
and fitness was monitored after 
spotting onto solid agar plates 
at 384 spot cell density for a to- 
tal of 16 passages. (B) Cultures 
were serially pinned at 1536 
colonies per plate and photo- 
graphed for 22 passages. 



22 passages 



96-well liquid senescence assay 

For the liquid senescence assay, 96-well plates that contained 60 
different yfgA estlA genotypes along with 36 independently generated 
his3A estlA "wild type control" strains were inoculated. Cultures in 16 
such replicate plates, each with independently generated double 
mutants, were grown in parallel (Figure SI shows the plate layout). 
In combination with estlA, we tested the effects of deleting 27 genes 
previously reported as affecting telomere length (Askree et al 2004; 
Gatbonton et al 2006; Shachar et al 2008), 15 affecting DNA repli- 
cation and repair (Collins et al 2007), and 18 other genes that were 
randomly chosen controls or of some prior interest (Table SI). The 
his3A estlA strains are suitable "controls" because, although HIS3 is 
disrupted with the KANMX marker, the cells express spHIS5 from the 
MA Ta- specific promoter and are therefore histidine prototrophs and 
expected to show no growth defects (Costanzo et al 2010). In total, 16 
yfgA estlA and 576 his3A estlA replicate cultures were examined in 
these liquid-based experiments. 

To follow growth in liquid, we used the procedure outlined in 
Figure 1A. Two -hundred microliter (200 julI) cultures were incubated 
for 22 hr at 30°C without shaking, and then shaken to resuspend the 
yeast cells. A sample from each culture was then diluted approxi- 



mately 1:70 into fresh media at each passage using a pin tool. To 
assess fitness at each passage, a separate sample from each culture 
was also diluted approximately 1:70 in water and spotted onto agar 
plates, at 384 spots per plate. Each of four 384-spot plates was in- 
cubated for 48 hr and photographed to determine cell density (Figure 
2). A 1:70 dilution was a pragmatic choice of dilution in the range 
used in low-throughput experiments. In low-throughput experiments 
where we carefully count cell number at each passage and dilute back 
to 5 x 10 5 cells per milliliter every day (Figure S3), the dilution ranges 
from 1:600 for healthy cultures to 1:7 when cells are deeply senescent. 

It is clear that the majority of his3A estlA strains, situated in the 
outer two rows/columns of each 384-spot plate, entered into and re- 
covered from senescence over 16 passages (Figure 2). The majority of 
his3A estlA strains grew well at passage 1, by passage 5 they grew very 
poorly, by passage 10 some had improved, and by passage 15 the 
majority had returned to a fitness similar to that at passage 1. At 
passage 15, some his3A estlA cultures did not grow at all. This pattern 
of growth is consistent with the hypothesis that the SGA procedure 
worked efficiently and that most, but not all, telomerase- deficient 
his3A estlA cultures entered into and recovered from telomere - 
initiated senescence. Recovery, in particular, was stochastic. To follow 
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Figure 2 Medium throughput 
liquid culture based senescence 
assay. 96-well cultures were 
spotted at 384 spots per agar 
plate, and plates photographed 
and analyzed after 48 hr 
growth. (A) Example images 
showing growth at passages 1 , 
5, 10, and 15. Four independent 
cultures of five genotypes are 
illustrated by the colored boxes. 
(B) Mean density profile (MDP) 
analysis for six genotypes for 16 
passages. Mean culture size and 
standard deviations are shown. 
All images can be browsed on 
our supplementary website (http:// 
research.ncl.ac.uk/colonyzer/ 
ChangSenescence/#i mages). 



the behavior of individual cultures, it is useful to browse through 
senescent time courses, with each time point shown as a stack, in files 
(.tiff) from our supplementary data website (http://research.ncl.ac.uk/ 
colonyzer/ ChangSenescence/ #images) . 

Because most gene deletions we chose to analyze affect DNA repair 
or telomere biology, most other strains (e.g., rad52A estlA and upf2A 
estlA) grew poorly in comparison with his3A estlA strains (Figure 
2A). RAD52 is required for recombination-dependent maintenance 
of telomeres in the absence of telomerase, and as expected, by the 
first passage the majority of rad52A estlA strains were growing poorly 
and never recovered (Mceachern and Haber 2006). The slight increase 
in fitness at the end of the time course for rad52A estlA strains (Figure 



2B) may reflect a £AD52-independent mechanism of survivor gener- 
ation that occurs in this and some other genetic backgrounds 
(Grandin and Charbonneau 2009; Lebel et al. 2009). Interestingly, 
the effect of upf2A, which affects nonsense-mediated RNA decay, has 
previously been reported to delay entry into senescence (Enomoto 
et al. 2004), whereas our data suggest that upf2A accelerates senes- 
cence in the estlA background and that upf2A estlA strains only 
weakly recover (Figure 2). 

chklA estlA strains behaved very similarly to his3A estlA strains, 
whereas, interestingly, riflA estlA and rif2A estlA strains recovered 
from senescence more rapidly than the bulk of other strains and the 
control his3A estlA strains. The effects of riflA and rif2A can be seen 



200 I H-Y. Chang et a/. 



.^5&G3' Genes | Genomes | Genetics 



from the high culture densities for these strains in passages 5 through 
10, relative to neighboring cultures (Figure 2A). This is of interest 
because Rifl and Rif2 have well-known roles in telomere biology, 
encoding proteins which bind via Rapl to telomeric DNA, and are 
known to affect telomere capping (Addinall et at 2011; Anbalagan 
et at 2011; Chang et at 2011). Interestingly, and consistent with our 
results, it was recently reported that rifl A and rif2A strains enter and 
recover from senescence more rapidly than other cell types (Chang 
et at 2011). The positions of all other genes tested by this liquid 
senescence assay are shown in Figure SI, and one complete series of 
photographs are shown in Figure S2A. To utilize data from all repli- 
cate cultures and to allow quantitative comparison between strains, we 
constructed mean density profiles (MDP) for all genotypes. MDPs are 
curves showing the senescence characteristics of each genotype and 
are constructed after estimating yeast culture cell densities for repeat 
cultures of different genotypes at each passage from photographs 
using the image analysis tool Colonyzer (Lawless et at 2010). Six 
examples of this type of analysis showing the MDP of each are shown 
in Figure 2B, and all 61 genotypes are shown in File SI. Analysis of 
576 independent his3A estlA cultures shows that at passage 1 they 
grew well, with an average fitness of around 4000 units (Figure 2b). As 
expected, strains became less fit with each passage, reaching a nadir of 
fitness at around passage 6, before recovering reasonable fitness by the 
end of the experiment at passage 16. The shape of this curve, com- 
bined with the fact that some cultures became sterile, is further evi- 
dence that the SGA procedures were efficient and that little 
contamination by other genotypes {e.g., diploids) or organisms (bac- 
teria) was observed. By passage 16, at the end of the experiment, the 
mean fitness of his3A estlA strains was nearly as high as that at 
passage 1, which was as expected because of the emergence of re- 
combination-dependent survivors; however, a large variance in cell 
densities was observed (Figure 2B). This variance at passage 16 was 
largely due to some cultures being essentially dead (see Figure 2A), 
presumably because in these particular cultures no senescent cells 
recovered. Sterile cultures were expected because the majority of 
cells do not recover from senescence and only a small volume of cell 
culture (<3 ul) was transferred each time. Variance also began to 
increase earlier in the experiment, most likely due the stochastic na- 
ture of entry into senescence. 

As expected from the examples shown in Figure 2 A, chklA estlA 
strains showed a similar MDP to his3A estlA strains, although in- 
terestingly, there was less variance between repeated chklA estlA cul- 
tures and less evidence of loss of vitality with time (Figure 2B, second 
panel from top). Quantitative analysis of 16 rad52A, upf2A, rifl A, and 
rif2A strains was consistent with the images. 

We used MPDs to visually classify each of the 60 gene deletions 
into different classes (File S3). These included those that entered 
senescence rapidly, including mrellA, xrs2A, radS2A, rad54A, radSSA, 
rad57A, rttlOlA, rttl07A, rttl09A, pol32A, upflA, upf2A, upf3A, dcclA, 
ctf8A, ctfl8A, sbc5A, sbc8A, sgslA, rad27A, mmslA, mms22A, cdc73A, 
rtflA, rfmlA, sumlA, kemlA, rrp8A, asflA, tellA, elglA, hmolA, edelA, 
ydlll8wA, rho4A, spt21A, tsalA, leal A, and slalA; a large number that 
behaved normally, including exolA, chklA, rad9A; and a small 
group of just three genes, est3A, rifl A, and rif2A, consisting of 
mutants that appeared to enter senescence early and recover early. 
Table SI summarizes the behavior of individual genes and whether 
the results we obtained were consistent with previously published 
studies. It is clear that our results largely reproduced results pre- 
viously reported from detailed, high-quality, low-throughput experi- 
ments giving us confidence in our data. However, there were some 
interesting surprises, such as the behavior of upflA, upf2A, and upf3A 



mutants that failed to recover from senescence here (Figure S5) but had 
previously been reported to delay senescence (Enomoto et at 2004). 

In summary, the 96-well senescence assay was capable of measuring 
entry of telomerase-deficient cells into senescence and recovery. 
However, this procedure was comparatively inefficient and not readily 
scalable to performing a genome-wide analysis of genes that affect entry 
into/exit from senescence. Therefore, we also tested a much more 
parallel and efficient solid pinning procedure, discussed next. 

1536-pin solid-agar senescence assay 

By pinning parallel colonies rather than growing liquid cultures, we 
increased throughput 36-fold. In pinning experiments, we grew a 
minimum of 8, sometimes 56 (for those on plate 15, which contained 
24 extra copies of a small number of gene deletions), and in the case of 
his3A, 144 independently generated double mutant cultures for a total 
of 22 passages. Comparing cell density estimates immediately after 
inoculation with estimates in stationary phase for robustly growing, 
1536-format, pinned cultures, we estimate that the maximum dilution 
achievable at each passage with this method is 16-fold. It is clear that, 
as cultures were passaged, the fraction of viable cultures decreased, 
suggesting cells entered into senescence (Figure 3A). In some cultures, 
culture size decreased to a minimum and then increased (with in- 
creasing passage), suggesting that these cultures recovered from senes- 
cence. Figure 3 A shows plate 15 that, in most cases, contains 24 
independently generated yfglA estlA double mutants and numerous 
his3A strains. Overlaid, colored boxes highlight the positions of his3A, 
upf2A, rad52A, and chklA strains, allowing a visual assessment of 
reproducibility between cultures. Figure S2 compares the growth of 
cells in the 1536 colony format with the lower throughput, 384-spot 
format across several passages. All senescent time course photographs, 
with each time point shown as a stack in files (.tiff), are in our 
supplementary data website (http://research.ncl.ac.uk/colonyzer/ 
ChangSenescence/#images). 

If strains had entered into and recovered from senescence, then 
telomere structure should alter, showing Type I- and Type II-type 
structures typical of cells using recombination-dependent ALT mech- 
anisms to maintain telomeres (Lundblad and Blackburn 1993; Teng 
and Zakian 1999). Importantly, Southern blots showed that the telo- 
mere structures of strains that had been passaged on solid plates had 
indeed altered, forming typical Type I or Type II survivors (Figure S4). 
In addition many of the DNA repair mutants known to be defective in 
recovery from senescence did not recover from senescence in these 
highly parallel experiments (e.g., rad52A estlA, Figure 3A). These two 
observations strongly suggest that most, if not all, estlA cells that were 
remaining after 22 passages had entered into and recovered from 
senescence and, therefore, that the high-throughput 1536 solid pin- 
ning approach is a useful method for following the senescence process. 

We generated MDPs of all 4299 genotypes over 22 passages, and 
plots of these profiles together with culture size standard deviations 
are available (File S2). It is interesting to compare the solid pinning 
data (File S2) with that for the 60 strains also analyzed in liquid 
culture, 6 of which are shown in Figure 3B, and the remainder of 
which are shown in File S2. In contrast to the pattern observed in 
liquid cultures, in pinning experiments evidence of a period of poor 
growth followed by recovery for many strains was less striking (e.g., 
MDPs for his3A estlA and chklA estlA strains, Figure 3B). However, it 
is notable that the variance in fitness increased with passage for each 
of these strains, more so for his3A estlA than for chklA estlA strains, 
and that his3A estlA strains showed an overall decline in fitness that 
plateaued around passage 15. We believe that the lack of clear decline 
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Figure 3 High-throughput, solid- 
culture-based senescence assay. 
1536 parallel colonies were 
transferred from agar plates, 
photographed, and analyzed. 
(A) Example images showing 
culture size at passages 1, 6, 
1 1 , and 22. Representative sets 
of four independent cultures of 
four genotypes are indicated by 
colored boxes. Plate number 1 5 
is presented, which has 24 in- 
dependently generated double 
mutant cultures per genotype, 
whereas plates 1 to 14 only 
have 4 independent cultures 
per strain. (B) Mean density 
profiles (MDP) for six genotypes 
for 22 passages. Mean culture 
size and standard deviations are 
shown. (C) Examples of rifIA 
and rif2A strain growth during 
passage. All images can be 
browsed on our supplementary 
website (http://research.ncl.ac.uk/ 
co I onyzer/Chang Senescence/ 
#images). 
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and recovery in the solid experiments is because of the low dilution at 
each passage (16-fold in solid vs. 70-fold in liquid). 

rad52A estIA strains entered senescence and failed to recover, as 
expected, but interestingly showed different growth patterns in solid 
vs. liquid plates. At passage 1, the fitness of rad52A estIA strains (1400 
units) was clearly less than the fitness of estIA chklA and estIA his3A 
strains (2000 units). This most likely is due the low fitness induced by 
deletion of RAD52 itself. With time, the fitness of rad52A estIA de- 
clined until, by passage 10, cultures were very small and did not really 
recover. The rate of entry of rad52A estIA cultures into senescence was 
much slower than that observed in liquid cultures, consistent with 



cultures being diluted less on solid agar compared with liquid cultures 
and, therefore, having less potential to divide. 

Curiously, the density profiles of upf2A estIA strains on solid pins 
seemed to be more similar to his3A estIA strains than to rad52A estIA 
strains. Conversely, in liquid cultures, upf2A estIA strains are more 
similar to rad52A estIA strains (compare the density profiles in Fig- 
ures 2B and 3B). We are unsure of the reason for this difference, but it 
seems to be reproducible as the effects of upflA, upf2A, and upf3A 
mirror each other in both solid and liquid cultures (Figure S5). Dis- 
rupting the nonsense-mediated RNA decay pathway affects many 
different transcripts and, hence, cellular processes, including those 
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Figure 4 Three senescence pro- 
file archetypes. Mean density 
profiles (MDP) were estimated 
for each of 4300 gene deletions 
combined with est7A. Pearson's 
correlation coefficient (r) was 
calculated for each mutation 
compared with his3A, rad52A, 
and rifIA mutations, represent- 
ing three distinct senescence 
curve archetypes. (A) A heat 
map showing the MDP of the 
mutations having the 1 4 highest 
correlations with that of his3A 
estIA. (B) All 3035 MDPs (pink) 
correlating strongly (r > 0.85) 
are shown along with the MDP 
of the archetype his3A est 7 A 
(red), and the mean curve of 
the 3036 strongly correlating 
MDPs (black). (C) A heat map 
showing the MDP of the muta- 
tions having the 20 highest cor- 
relations with that of rad52A 
estIA. To construct this group, 
we only included strains whose 
culture size at passage 2 was 
less than C s = 1680, where 
C s = C ave - 0.5*C SD , C ave is 
the average of all double mu- 
tant culture sizes at passage 2, 
C SD is the culture size standard 
deviation at that passage. (D) 
All 208 MDPs (pink) correlating 
strongly (r > 0.5, culture size < 
C s ) with the MDP of rad52A 
estIA (red), together with the 
average of the 209 strongly cor- 
relating MDPs (black). (E) A heat 
map showing the MDP of the 
mutations having the 14 highest 
correlations with that of rifIA 
estIA. (F) Twenty gene disrup- 
tions whose MDPs (pink) corre- 
late best with the MDP of rifIA 
est7A (red), together with the 
average of the 21 most strongly 
correlating MDPs (black). (G) 
Scale for heat maps. 



relevant to telomere biology (Addinall et ol. 2011; Enomoto et al 
2004). Thus the explanation for the differences between solid and 
liquid cultures may be as trivial as the upfA strains forming cultures 
more easily from dense inoculum after pinning on solid agar than 
from the dilute inoculum used in liquid cultures. Finally rifIA estIA 
and rif2A estIA mutants essentially behaved as expected given the 
liquid culture results and, if anything, improved fitness with time 
(Figure 3B). 

Classifying senescence curves: three archetypes 

To help classify the behavior of the thousands of different genotypes, 
we generated MDPs for each of the ~4300 gene deletion strains (File 



S4). We classified strains according to the similarity of their MDP to 
that of three archetypal gene disruptions: his3A estIA, a neutral con- 
trol strain (Figure 4A, B and File S5); rad52A estIA, affecting a typical 
gene required for recovery from senescence (Figure 4C, D, File S6, and 
File S7); and rifIA estIA, showing accelerated entry into and exit from 
senescence (Figure 4E, F and File S8). Similarity between MDPs was 
quantified using Pearson's correlation coefficient (r). 

Normal senescence: his3A estIA archetype 

As we expected, the vast majority (more than 3000) of the double 
mutants can be classified in the category of "normal senescence" 
(Figure 4A), defined as mutants whose MDP correlates well with that 
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Figure 5 Properties of gene deletions that cause similar 
effects to rad52A in an est7A background. (A) Venn di- 
agram showing gene deletions reported to result in sig- 
nificantly shorter (red) or longer (green) telomeres 
(Askree et al. 2004; Gatbonton et a/. 2006; Shachar 
et al. 2008) and those we classified as similar to the 
rad52A estIA archetype (black). (B) Gene ontology (pro- 
cess) analysis demonstrating the P values for signifi- 
cance of over representation of terms (both GO terms 
and classifications from other screens) over represented 
in the class similar to the rad52A estIA archetype. 



of his3A estIA strains (r > 0.85). There was a general trend of lower 
density across all passages for these strains (Figure 4B); at late pas- 
sages, this is essentially due to mean density being estimated from 
a mixture of strongly growing cultures and sterile or poorly growing 
cultures {e.g., Figure 3A, P22 Supp Fig S2). In this class, the mean 
density plateaued around passage 18. In normal senescence strains, on 
average 50% of cultures recovered growth by passage 19. Interestingly, 
the gene deletions whose MDP have the strongest anticorrelation with 
his3A estIA strains include many genes affecting telomere mainte- 
nance, DNA damage response, and cell-cycle control (bottom of File 
S5). In fact, the five strains whose MDP is least well correlated with 
those of his3A estIA (rif2A estIA, rad27A estIA, riflA estIA, elglA 
estIA, and est3A estIA) are the strongest examples of the riflA estIA 
archetype, which we define below. 



Accelerated entry into and exit from senescence: riflA 
estIA archetype 

Using a strict correlation cut off (r > 0.85) leaves only three members 
of the riflA estIA class; however, we can see in Figure 4F that the 20 
most highly correlated genes display a distinctive MDP. Deletions in 
this class appear not only to speed up senescence but also to accelerate 
the appearance of survivors. For strains in this class, on average 81% 
of cultures recovered growth by passage 19. Interestingly, EST3, RIF1, 
RIF2, and TEL1 are genes that are well known to affect telomere 
maintenance and Chang et al. (2011) reported that riflA and rif2A 
mutants accelerated entry into and recovery from senescence in riflA 
est2A and rif2A est2A mutants. The effects of est 3 A on estIA mu- 
tants that we have observed are not consistent with previous reports 
(Lendvay et al. 1996). However, we have confirmed this phenotype in 
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Table 1 GO term analysis of RAD52-\\ke genes in this study 



Macromolecular complex terms 



Genes annotated to the term 



Ribonucleoprotein complex 
Ribosome 

Chromatin remodeling complex 

Ubiquitin ligase complex 

Preribosome 

SNARE complex 

Elongator holoenzyme complex 

Histone deacetylase complex 

Microtubule 

Chromatin assembly complex 
Histone acetyltransferase complex 
Ctf18 RFC-like complex 
MRX complex 
DUBm complex 
Rhp55-Rhp57 complex 
GSE complex 
Rpd3L complex 

Replication fork protection complex 
HOPS complex 
SetlC/COMPASS complex 
Mediator complex 
Vesicle coat 



CKB1, EAP1,MRPL36, MTC1, NOP12, RPL21B, RPL35B, RPS11B, 

RPSUA, RPS18B, RPS4A, TEF4, TMA19 
MRPL36, MTC1, RPL21B, RPL35B, RPS11B, RPSUA, RPS18B, RPS4A, 

TEF4, TMA19 
ARP6, HDA1, NHP10, SDS3, SIN3, SWD1, SWD3 
DIA2, MMS1, MMS22, NPL4, SLX5, SLX8, UBX2 
CKB1, NOP12, RPL35B, RPS11B, RPS16A, RPS4A 
GOS1, TLG2, VAM3, VAM7 
ELP2, ELP3, ELP4, ELP6 
HDA1, SDS3, SIN3 
CIK1, KAR3, SHE1 
ASF1, CAC2, HTZ1 
SGF11, SGF29, SGF73 
CTF18, CTF8, DCC1 
MRE11, RAD50, XRS2 
SGF11, SGF73 
RAD55, RAD57 
GTR1, MEH1 
SDS3, SIN3 
CTF4, MRC1 
VAM6, VPS41 
SWD1, SWD3 
MED 7, SOH1 
LST4, SEC28 



RAD52-\\ke genes in this study were analyzed using SGD gene ontology slim mapper (http://www.yeastgenome.org/cgi-bin/GO/goSlimMapper.pl). "Macromolecular 
complex terms: Component" GO set was used to determine the list of protein complexes. 



careful low throughput experiments in the different W303 genetic 
background (Figure S3). The tell A estlA results from previous studies 
are somewhat conflicting (Abdallah et al 2009; Chan et al 2001; 
Enomoto et al 2002; Ritchie and Petes 2000, Table SI). Some previous 
data agree with our results whereas other data do not; it is not clear 
whether the differences are due to strain background or some other 
factor, such as experimental design. Other genes in this class, notably 
ELG1 encoding a replication factor C (RFC) -type protein and with 
known telomere-related phenotypes (Smolikov et al 2004), warrant 
further study. 

No recovery from senescence: rad52A estlA archetype 

There were 208 double mutants (File S7) having an MDP that corre- 
lated well with that of rad52A estlA (r > 0.5, culture size at passage 
2 < 1680, Figure 4C, D). For strains in this class, on average 29% of 
cultures recovered growth by passage 19. This subset is interesting 
because it contains many genes required, in principal, for recovery 
from telomere-initiated senescence, and orthologs of such genes 
could play critical roles during carcinogenesis, particularly in ALT 
tumors. Consistent with the role of this gene class at telomeres, 
there was a strong enrichment of gene deletions previously found 
to result in short telomeres but not in long telomeres (Askree et al 
2004; Gatbonton et al 2006; Shachar et al 2008; Figure 5A). This 
enrichment is perhaps not surprising because cells with gene dele- 
tions inducing short telomeres are likely to enter senescence faster 
than those with gene deletions that result in normal length or long 
telomeres. 

To further characterize this class, we subjected the genes with 
MDP similar to the rad52A estlA archetype to a "gene ontology (GO) 
yeast process" analysis. Interestingly, we found significant enrichment 
for genes involved in response to stress, chromosome organization, 
DNA metabolic processes, cellular membrane organization, vesicle- 
mediated transport, transcription DNA dependence, cell cycle, and 



protein modification (Figure 5B). Our results suggest that these genes 
share many properties and may play an important role in replicative 
senescence and might do so by affecting several cell processes. More- 
over, when subjecting these genes to a GO analysis, we identified 
many components of known and relevant complexes, for instance, 
the MRX complex and the Ctfl8 RFC-like complex (Table 1). The 
MRX complex has many known roles at telomeres, and when 
deletions are combined with deletions in the Ku complex, the 
double mutants enter senescence like telomerase-deficient strains 
(Dubois et al 2002; Maringele and Lydall 2004a). However, we are 
unaware of reports of the MRX complex inhibiting senescence. The 
Ctfl8 complex does have roles at telomeres (Addinall et al 2011; 
Hiraga et al 2006; Khair et al 2010); however, as far as we are 
aware, this is the first report showing that deletion of components 
of this complex accelerates senescence. The mechanisms by which 
these complexes affect replicative senescence need to be clarified in 
future studies. 

Genes similar to the rad52A archetype have negative 
genetic interactions with yku70A but not cdc13-1 

The high-throughput senesence assays we report here are comple- 
mentary to other high-throughput studies aimed at understanding 
telomere or chromosome biology. For example, Addinall et al (2011) 
performed systematic suppressor/enhancer analyses on yeast strains 
containing temperature-sensitive telomere capping defects caused by 
cdcl3-l or yku70A mutations. All gene deletions were classified as 
suppressing, enhancing, or showing no strong interaction with each 
defect, and interaction strengths were quantified. The two datasets 
were compared in a single plot, and it is interesting to overlay gene 
deletions that behave like rad52A estlA in high- throughput senescence 
assays on this plot (Figure 6A). 

Interestingly, most of the gene deletions that accelerated senes- 
cence also enhanced yku70A temperature-sensitive growth defects but 



:~5iG3'Genes | Genomes | Genetics 



Volume 1 August 2011 I Regulation of Senescence I 205 



enhancement 



suppression 



LST4 

RPN4 



MEH1 



0) 

c 

o 




. . . AR02 i 



EL*M1 



VAM10 - • 

R lftpL024W # 

YPT7 «fo»„ "ffi? S ^jf xl2 SN ' M 

■EST ** u SLX8 Hftw*"^,. I s * . 




6 



- 0.5 0.0 
Genetic interaction score: yku70A. 37.5° 

• Gene deletions that behave similarily to rad52A in this study 

• Gene deletions that behave similarily to rad52A in this study, also have short telomere phenotype 
Gene deletions that behave similarily to rad52A in this study, also have long telomere phenotype 



B ft 



Deletions that behave similarily 
to rad52A in this study (209 genes) 



All genes (4294) 



76% I 24% 



57% 



43% 



44% I 56% 



50% 



50% 



Figure 6 Gene deletions that 
cause similar effects to rad52A 
in an est7A background are 
enriched for deletions that neg- 
atively interact with yku70A but 
not cdc13-1. (A) Gene deletions 
identified in this study as accel- 
erating senescence overlaid on 
data showing genes that sup- 
press or enhance temperature- 
sensitive growth defects caused 
by cdc13-1 or yku70A. The data 
points are from Figure 4 in 
Addinall et a/. 201 1 . Genes are 
highlighted in blue if their mean 
density profile (MDP) is classed 
as similar to the rad52A estIA 
archetype, in red if they addi- 
tionally have a short telomere 
phenotype, or in green if they 
additionally have a long telo- 
mere phenotype as defined be- 
fore (Askree et a/. 2004; 
Gatbonton et a/. 2006; Shachar 
et a/. 2008) (B) Left: Fraction of 
strains similar to the rad52A 
estIA archetype found in indi- 
cated regions of the genetic in- 
teraction profile presented in 
panel A. Right: Fraction of total 
gene deletions tested found in 
indicated regions of the genetic 
interaction profile presented in 
panel A. 
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showed no significant bias toward either positive or negative effects on 
the fitness of cdcl3-l mutants. This is shown by the clear left/right but 
lack of top/bottom bias of such strains in Figure 6. The yku70A 
mutants have short telomeres, similar to those of telomerase-deficient, 
pre-senescent cells, whereas cdcl3-l mutants have telomeres of nor- 
mal length or slightly longer than normal (Zubko and Lydall 2006). 
Therefore, a simple explanation for this bias is that the gene deletions 
highlighted on the left of Figure 6A contribute to the ability of cells 
with short telomeres to continue through the cell cycle. These data 
also suggest that yku70A is perhaps a better model for telomere-ini- 
tiated senescence than cdcl3-l. 



Genes required for efficient replication 
fork progression 

Global genetic analyses in budding yeast Saccharomyces cerevisiae 
have identified a pathway, including ASF1, RTT109, MMSi, 
MMS22, RTT101, and RTT107, instrumental to cell survival after 
replisome stalling (Collins et al 2007; Duro et al. 2008; Pan et al. 
2006; Zaidi et al 2008). Interestingly, it has been proposed that 
RAD52 and MMSI are essential for pre-senescence growth (Abdallah 
et al 2009). To further investigate this pathway in telomerase-deficient 
cells, we examined the behavior of these strains in both our high- 
throughput and low- throughput screens. We found that not only 
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RAD52 and MMS1 but also ASF1, RTT109, MMS22, RTT101, and 
RTT107 showed an accelerated senescence phenotype when deleted 
in both screens (Figure 2, Figure 3, Figure S6). These results strongly 
suggest that replication fork progression plays an important role in 
replicative senescence. However, our data also suggest that there are 
many additional pathways/complexes affecting replicative senescence 
(Figure 5, Table 1). It will be interesting to investigate the roles of 
these pathways/complexes further. 

DISCUSSION 

We report the first genome-wide analysis of the effects of gene 
deletions on entry into and exit from telomere-initiated senescence in 
S. cerevisiae. We have successfully carried out a medium-throughput, 
liquid-based senescence screen and a high-throughput, solid pin-based 
screen. Although the liquid-based screen is a valuable source of data 
describing the behavior of 60 gene deletions hand-picked for the 
relevance of their effect on entry into and recovery from senescence, 
we chose to perform a genome-wide analysis of 4299 genes by using 
a solid pinning approach because it was much more efficient (36 times 
the throughput). Inevitably, the high-throughput pinning approach 
suffers from a lack of resolution compared with the liquid-based 
approach. Nevertheless, it contains a great deal of valuable quantita- 
tive information. We have constructed a high-quality dataset demon- 
strating senescence profiles for all 4299 gene deletions in the knockout 
library combined with a telomerase defect that will complement other 
low- and high-throughput analyses of telomere biology and the DNA 
damage response. 

We chose to classify double deletion mutations based on the 
similarity of their MDPs to those of archetypal mutations: his3A estlA 
(normal), rad52A estlA (accelerated entry into senescence and failure 
to recover), and riflA estlA (accelerated entry into senescence and 
recovery). However, other clustering or analysis approaches may also 
prove informative. For example, we also performed unsupervised hi- 
erarchical quality threshold (QT) clustering to identify 23 different 
MDP classes (Figure S8). The MDPs which we constructed (e.g., File 
S2 and File S4), together with the images and raw quantifications from 
which they are derived, are published to allow others to analyze our 
dataset using different approaches, for example, by searching for 
strains whose deletions result in MDPs correlating strongly with those 
of archetypes of specific interest, which could be different to the three 
presented here. 

The rad52A estlA archetype strain and mutants in this class show 
poor growth when generated by SGA immediately after germination. 
Therefore we tested whether additional information is gained from 
continuing to passage such poor-growing strains. To address this 
issue, we chose to compare our rat/52A-like cluster with genes we 
identified as growing poorly with estlA at passage 0 in our experi- 
ments, and with analogous datasets from large-scale, genome-wide 
synthetic lethal analysis experiments (Costanzo et al. 2010). Figure 
S7A shows that there is overlap among these three datasets but that 
they are clearly not identical. The group of strains with MDP similar 
to that of the rad52A estlA archetype is distinct from that which we 
identified as synthetic lethal interactions immediately after SGA (File 
SlOandFileSll). 

To compare the relative merits of each dataset, the member strains 
of each are highlighted in Figure S7C-F on plots similar to Figure 6. 
Interestingly, the left right/bias we observed for the strains in the 
rad52A estlA class (Figure 6 and S7C) is not as pronounced in the 
other datasets. This strongly suggests that gene deletions we identified 
in the rad52A estlA class, based on their behavior with repeated 
passage, have different properties than those growing poorly at early 



passages identified by us or in the Costanzo et al. (2010) dataset. 
Therefore we conclude that, at the very least, our high-throughput 
senescence assays complement the existing datasets, identifying gene 
deletions with notably different properties. 

Our dataset is a rich resource for those working in telomere 
biology in budding yeast and other eukaryotic organisms. All our data 
are freely available, and there are numerous provocative and inter- 
esting observations worthy of further investigation. For example, the 
genes behaving like RAD52 or RIF1 show very interesting effects on 
the senescence process in budding yeast. Although some genes in 
the RAD52 class, such as SLX5, SLX8, and MMS1, were previously 
reported to affect the processes of recovery from senescence, we have 
identified 200 or more genes in this class. Our data also shed new light 
on the functions of proteins with known roles in telomere biology, 
such as Elgl or members of the Ctfl8 complex. 

We imagine that iterations of our genome- wide approach to under- 
stand senescence, perhaps utilizing genome-wide liquid passage ap- 
proaches, like Figure 2, or "bar code-based" approaches (Hillenmeyer 
et al. 2008), may extend our systematic analysis of gene deletions 
affecting telomere-driven senescence. It may also be interesting to 
perform genome-wide screens analyzing the effects of deleting other 
protein or RNA subunits of telomerase, as we have shown that de- 
leting the EST 3 gene affects the senescence profile of estlA mutants in 
a manner analogous to riflA mutants. It may also be interesting to 
extend this approach to S. pombe, where telomere biology is signifi- 
cantly different from S. cerevisiae and genome- wide knockout libraries 
and SGA approaches are now available (Dixon et al. 2008; Roguev 
et al. 2007). It will be more technically challenging to perform anal- 
ogous studies in mammalian cells. 
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